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Abstract: We suggest a thick braneworld model in the squared curvature gravity theory. 
Despite the appearance of higher order derivatives, the localization of gravity and various 
bulk matter fields is shown to be possible. The existence of the normalizable gravitational 
zero mode indicates that our four-dimensional gravity is reproduced. In order to localize 
the chiral fermions on the brane, two types of coupling between the fermions and the brane 
forming scalar is introduced. The first coupling leads us to a Schrodinger equation with a 
volcano potential, and the other a Poschl- Teller potential. In both cases, the zero mode 
exists only for the left-hand fermions. Several massive KK states of the fermions can be 
trapped on the brane, either as resonant states or as bound states. 

Keywords: Field Theories in Higher Dimensions, Large Extra Dimensions 



1 The corresponding author. 



Contents 



1 Introduction 



1 



2 Model and solution 



3 



3 Tensor perturbations and the localization of gravity 



4 



4 Trapping of matter fields 

4.1 Massive fermions as resonant states 

4.2 Massive fermions as bound states 



9 

11 
11 



5 Conclusions and prospect 



13 



1 Introduction 

The suggestion that our world might be a four-dimensional domain wall was proposed early 
in the 1980s [1, 2], but it wasn't until the localization of gravity on either thin [3, 4] or 
thick [5-9] domain wall was shown to be possible, that the braneworld models attracted 
renew attentions. More and more studies indicate that the braneworld models might be 
effective in solving some long existed problems, such as the hierarchy problem and the 
cosmological constant problem [10-12]. In this paper, we focus mainly on the thick domain 
wall case. A comprehensive review about thick brane models can be found in [13]. 

Some questions needed be addressed in the thick braneworld models: the stability of 
the solution and the possibility of localizing matter fields as well as gravity on the brane. 
In the frame of Einstein's general relativity, there are a lot of works which devoted to 
answer these questions, see [5-7, 9] for the localization of gravity, and [14-25] for trapping 
various kinds of bulk matter fields on both thin and thick branes. For recent developments 
see [26-34]. 

However, due to the fact that general relativity is not renormalizable, we are suggested 
to consider the effects of higher order curvature terms [35]. Such terms also appear in 
the low-energy effective action of string theory [36]. To evade the appearance of spin-2 
ghosts, and the Ostrogradski instability [37], the higher order curvature terms usually are 
introduced as the Gauss-Bonnet term or an arbitrary function of the curvature, namely, 
f(R)- Both cases were applied to consider various of issues in cosmology and higher energy 
physics, for details, see [38-40] and references therein. 

As to the braneworld scenarios, it is shown in many works that the introduce of the 
Gauss-Bonnet term usually imposes no impact on the localization of gravity [41-50] . This 
result might change in f(R) gravity, however, because now some higher order derivatives 
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evolve in. The appearance of higher order derivatives usually renders the constructing of 
a braneworld models very hard. 

For pure gravitational systems which contains no any matter fields, some thin braneworld 
models have been constructed in a lower order frame [51-54] by using the Barrow-Cotsakis 
theorem [55]. This theorem states that the fourth order f(R) gravity (higher order frame) 
is conformally equivalent to a second-order gravity theory (lower order frame). Unfortu- 
nately, for thick domain wall models, this method would lead to ambiguous [56] and thus 
is not a good choice. 

Some thick brane solutions were found directly in the higher order frame [57, 58]. 
In [57], with a background scalar field, the authors found some analytical thick brane 
solutions in both constant and variant curvature cases. While in [58], by analyzing the 
existence of the fix points, the authors found some numerical thick /(i?)-brane solutions 
with pure gravity. The trapping of complex scalar field on the brane solutions was shown 
to be possible in [58]. Both papers considered the case where f(R) oc R n . 

However, the solutions given in [57, 58] are not perfect for the following reasons: 

• The solution in M5 space given in [57] contains singularity, because the warp factor 
is not smooth at the brane. Meanwhile, even though the solutions found in both dS§ 
and AdS§ spaces are stable under the tensor perturbations, only the solution in dS§ 
space supports normalizable zero gravitational mode, see [59]. Therefore, the gravity 
can be localized only on the brane in dS§ space. However, the constant curvature 
space is a rather special case. 

• For the case of non-constant curvature, some solutions were suggested in [57, 58]. 
However, these solutions are not very good: the solution in [57] contains a singular 
point at where the curvature diverges; while the solution in [58] is given numerically, 
and we only know its behavior in infinity. 

• Besides, the model in [58] is purely gravitational, namely, no matter field were intro- 
duced. However, to localize fermions, one usually needs to couple the fermion with a 
scalar field. 

For the reasons above, we aim at constructing a domain wall brane which has the 
following properties: 

• The brane solution is found in a system which contains higher order curvature terms 
and especially a system contains fourth-order derivatives. 

• The brane is smooth, stable and is formed by a background scalar field. 

• The gravity can be localized on such brane. 

• By introducing the Yukawa coupling between the fermion and the background scalar 
field, the fermions can be trapped on the brane. 
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As a toy model, let us focus on the squared curvature gravity, i.e., the higher curvature 
correction is given by f(R) oc R 2 . For this special gravity, a thin braneworld has been 
constructed in [51] by using the Barrow-Cotsakis theorem. 

In the next section, we give a thick domain wall solution in the squared curvature 
gravity. In section 3, we discuss the localization of four-dimensional gravity. In section 4, 
the trapping of bulk matter fields with various spins is analyzed, but we mainly focus on 
the half spin fermions. Our conclusions are listed in the last section. 



2 Model and solution 

We start with the action 



1 1 



S = I d^x^g ( ^f(R) - -d M ct>d M ct> - V{<j>) ) , (2.1) 



2k% ' v ' 2 

where k| = 87rG5 is the gravitational coupling constant with G5 the five-dimensional 
Newtonian constant. We have introduced a background scalar field <f>, of which the self- 
interacting is described by the potential V(4>). We prefer to use Roman letters to denote 
the bulk coordinates, i.e., M,N = 0,1,2,3,4 and the Greek letters represent the brane 
coordinates, i.e. /i, v = 0, 1, 2, 3. 

Our discussions will be limited on the static flat braneworld scenario, for which the 
metric is given by 

ds 2 = e 2A ^r]^dxi*dx v + dy 2 , (2.2) 

with y = x 4 the extra dimension. Meanwhile, the scalar field, (f> = 4>(y), is independent of 
the brane coordinates. For system (2.1)-(2.2), the Einstein equations read 

f(R) + 2f R {4A' 2 + A") - 6f' R A' - 2/£ = k 2 ^' 2 + 2V), (2.3) 

and 

- 8f R {A" + A' 2 ) + 8f R A> - f(R) = 4(0> 2 - 2V), (2.4) 

where the primes represent derivatives with respect to the coordinate y and f R = df(R) /dR. 
The scalar field is described by the following equation: 

8V 

4AV + </,"- — = 0. (2.5) 

We consider a toy model, in which f(R) = R + ^R 2 . Then we would have f R = 
l+2 7 i?, f' R = 2-fR', f R = 2-fR", R! = -(8A"'+40A'A"), R" = -(8A""+40A" 2 +40A'A"'). 
When 7 = 0, we come back to the case of general relativity. Obviously, if 7 7^ 0, the Einstein 
equations are fourth-order differential equations. In order to find a domain wall solution, 
let us consider the following (j) 4 model: 

V(cj ) ) = \( 5 \^-v 2 ) 2 + A 5 , (2.6) 
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where A^ 5 ** > is the self coupling constant of the scalar field, and A5 is another constant. 
Obviously, at <j> = ±v the scalar potential V(4>) takes its minimum value, i.e., V(4> = = 
A 5 - 

It is easy to proof that eq. (2.5) supports the following solution: 



4>(y) = utanh 



e A ^ = cosh" 

However, the Einstein equations (2.3)-(2.4) strongly constrain the possible values of the 
parameters. Specifically, the Einstein equations are satisfied only when the parameters 
take the following values: 

A<»_-M, A 5 = -iI^. (2.9) 

784 7 ' y 20k 2 ' 6728 7 k§ V ; 

Defining a new parameter by 




(2.10) 



we can rewrite (2.7)-(2.8) in a more simpler form: 

(f>( y ) = vtanh(ky), e A(y) = cosh(fcj/) -1 . (2.11) 

For this reason, in the discussions below, we always use the parameter k rather than 7, to 
reflect the influence from the squared curvature term. 

As shown in figure 1, the scalar field is a kink which satisfies (f)(0) = 0, </>(=Loo) = ±v. 
V(4>) takes its minimum at (f) = ±v. We also see how the energy density p = Too = 
e 2A (^0' 2 + V ((/>)) distributes along the fifth dimension. Obviously, p peaks at y = 0, 
where the brane locates at. It is not difficult to verify that the geometry of the space- 
time becomes anti-de Sitter at y = ±00, where the bulk curvature R = —20k 2 and the 
corresponding cosmological constant is nothing but the parameter A5. The value of A5 can 
also be obtained from eq. (2.3) by taking the limit y — > 00. 

3 Tensor perturbations and the localization of gravity 

To discuss the localization of massless four-dimensional gravity, we need to consider the 
tensor perturbations: 

ds 2 = e 2A ^ (^ + h^dx^dx" + dy 2 , (3.1) 

where h^ u = h fiu (x p , y) depend on all the coordinates and satisfy the transverse-traceless 
condition 

»rv = = dj? v . (3.2) 
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The perturbation of the scalar field is denoted as 8<j) = (p(x^, y). 

In [59], we have demonstrated that the scalar perturbation <p decouples from the ten- 
sor perturbations h^ iu under the transverse-traceless condition (3.2), and the perturbed 
Einstein equations are given by 



d 2 y +^d y + e- 2A U^)h, u 



}r 



^8 h 



or in a more simpler form 



flip) L 
JR 



(3.3) 



(3.4) 



rfd^dv and CK 5 ) = g MN V mV n are the four-dimensional and five-dimensional 



Here CK 4 ) = if v 
d'Alembert operators, respectively. 

Obviously, for R = const, or f{R) = R, which simply gives the theory of general 
relativity, eq. (3.4) reduces to the five-dimensional Klein-Gorden equation for the mass- 
less spin-2 gravitons, and the reproduction of brane gravity has been discussed in [5-7]. 
However, for an arbitrary form of f(R) with variant R, the issue of localization of gravity 
should be given a further consideration. 

Following the method given in [5-7], we introduce a coordinate transformation dz = 
e~ A ( y ^dy, under which the metric (2.2) is conformally flat: 



aV 



^(rj^dx^dx" + dz 2 ). 



(3.5) 



It is very easy to proof that in the new coordinates system, y = arcsinh(fez)/fe and therefore 
the warp factor takes the simple form: 



3 2 ^) = (l + k 2 Z 2 )- 1 



Then the perturbed equation (3.4) reads 

d 2 



+ 



3^ + %^U + n( 4 ) 

a Jr 



0. 



(3.6) 



(3.7) 



Here we have denoted a(z) = e A ^ to write the final result looks more symmetrical. 



By doing the decomposition h^ v {x p ,z) = (a 3 ^ 2 f R )e^ v {x p )tp{z), with e^ u {x p ) satis- 

= d^u, we obtain a Schrodinger like 



fying the transverse-traceless condition rj 
equation for i/j(z): 

[d 2 - W(z)] </>(*) 
with the potential given by [59] 



-m 2 ^(z) 



W(z) 



3 a' 2 3 a" 3 a'f' 



4 o 2 + 2 a + 2 af R 



4/1 2f R 



(3.8) 



(3.9) 



One can also factorize the Schrodinger equation (3.8) as 



,'3d z a ld z fR 
l 2- + 2^f 



3 <9 2 a l d z f R 
2 a 2 f R 



iP( z ) = - m 2 ip(z) (3.10) 
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to conclude that there is no gravitational mode with m 2 < and therefore the solution in 
(2.7)-(2.9) is stable. For eq. (3.8), the zero mode (for which m? = 0) reads 

^(z) = N a 3 / 2 (z)f 1 R /2 (z), (3.11) 

where No represents the normalization constant. 
By plugging (3.6) into (3.9), one obtains 

15k 2 (-14 + 37k 2 z 2 + 28/eV + 4k 6 z 6 ) 

W(z) = * - '-, 3.12 

V ' i(5 + 7k 2 z 2 + 2k 4 z 4 ) 2 V ; 

the famous volcano potential [3, 5-7]. This potential supports a normalizable zero mode 
(see figure 3) 



<k V5 + 2k 2 



2 



^°)(z) = JOL v °^ * (3.13) 

l/ 8(1 + ^2)5/4' V ) 



as well as a series of continuous massive KK modes. 

Starting from m 2 > 0, the continuum of KK modes might lead to a correction to the 
Newtonian potential on the brane. As depicted in figure 2, W{z) ~ 4^ as \z\ 3> 1. It is well 
known (see, for example [6, 9]) that if the potential W(z) ~ a(a + l)/z 2 as \z\ 3> 1, then 
the KK modes on the brane take the form ^("^(O) ~ ?n a_1 , and for two massive objects 
at a distance r from each other, the correction for the Newtonian potential is AU oc l/r 2a . 
For our model a = 3/2, thus \^ m ^ (0)| 2 ~ m for small masses. Thus the correction to the 
Newtonian potential is AU oc 1/r 3 , which is same to the one given by the RS model [3]. 

Note that for KK modes with < m 2 < W m3iX , there might exist some resonant modes 
due to the quantum tunneling effect [5, 6]. Here W umx is the maximum of the potential 
W(z). To discuss this issue, let us redefine the variables as 

z = kz, m 2 = m 2 /k 2 . (3.14) 

As a consequence, one can rewrite equation (3.8) as 



d 



15 (-14 + 37 z 2 + 28z 4 + 4z 6 ) 



4 (5 + 7z 2 + 2z 4 ) 



ip(z) = -fh 2 ij(z), (3.15) 



which is parameter independent. This equation is hardly to be solved analytically, but we 
can apply some numerical methods to analyze the resonant structure. 

To start with, let us impose some initial conditions to simplify the numerical calcula- 
tion [30, 60-62]: 

VWen = Cl, ip' cvcn = 0; Vodd = 0, V>odd = c 2- (3-16) 

Here we have denoted VWem ^odd as the even and odd parity modes of i^{z), respectively. 
The prime represents the derivative with respect to z. In principle, c\ and C2 are arbitrary 
constants, here we simply take c\ = Oi = 1. The massive modes ip(z) here are analogous to 
the wave function in the quantum mechanics. Therefore we can interpret \ip(z)\ 2 dz (after 
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Figure 2. Plots of W(z) and z 2 W(z). 



normalizing ip(z)), as the probability for finding the massive KK modes within (z, z + dz). 
However, the massive modes are non-normalizable because they oscillate strongly when 
z> 1. Following the method given in [30, 60-62], we define the function 

p ^ 2) = M^Wi <3 - i7) 

as the relative probability for finding a massive KK mode in a narrow range —z\> < z < Zf, 
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Figure 3. The shapes of the gravitational zero mode ip ( z )- 

(as compared to a wider interval — z max < z < z max ) with mass square m 2 . Here 2zb is 
about the width of the thick brane and z max = lOzj,. For an eigenvalue m 2 3> W max , the 
corresponding KK mode ip{z) can be approximately identified as a plane wave, for which 
the probability Po,{fh 2 ) tends to 0.1 . 

According to the numerical calculation of Po,{fh 2 ), see figure 4, there exists no resonant 
state for the Schrodinger equation (3.15). Obviously, as the eigenvalue fh 2 approaches to 
the maximum of the scaled potential W max (z) = 1.1086, the relative probability Pq tends 
to 0.1, as expected. 
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Figure 4. The relative probability Pq for various values of m 2 
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The analysis above indicates that the solution (2.7)-(2.9) is stable under tensor pertur- 
bations, and the four-dimensional massless gravity can be localized on the brane. Besides, 
there exists no resonant gravitational KK modes. Thus the massive gravitons can not be 
quasi-localized on the brane. The continuous KK spectrum contributes a correction to the 
brane Newtonian potential, i.e., AU oc r~ 3 . Now let us move on to the issue of trapping 
matter fields. 

4 Trapping of matter fields 

An interesting but also vitally important question in brane world scenarios is how to localize 
various bulk matter fields on the brane by a natural mechanism. In many braneworld 
models the massless scalar fields can be trapped on the branes [26, 33, 63]. Usually, Spin-1 
Abelian vector fields can not be localized on five-dimensional flat branes. But, it can be 
localized on the RS brane in higher dimensional case [64] or on the thick dS brane [65] and 
Weyl thick brane [26, 28]. 

As far as our solution (3.6) is concerned, the scalar fields can be localized on the brane, 
while the vector fields can not. This conclusion can be achieved by following the discussions 
given in [28, 63]. 

In addition to the scalar and vector fields, the localization of Kalb-Ramond fields and 
gauge fields on single- or multi-branes were also investigated [31, 66-72]. But we are not 
going to discuss these issues. In stead, we would like to focus on the localization of the 
spin-1/2 fermions. 

According to the works in the past, fermions can not be localized on branes in five and 
six dimensions without introducing the scalar-fermion coupling [20, 24, 29, 30, 62, 73-77]. 
This is the reason why we insist on introducing a background scalar field. 

We proceed with the action 

= J d 4 xdyV^hr M (d M + UM)y-vyF(<t>)y\ (4-1) 

where F((f>) determines the coupling type and rj tells the strength. As in refs. [27, 28, 32], 
we denote T M = e~ A (-~f^, 7 5 ) as the five-dimensional gamma matrices with 7^ and 7 s the 
usual four-dimensional Dirac gamma matrices, ojm are the spin connections, of which the 
non-vanished components read = ^(dzA)"/^^. With the action in (4.1), we obtain the 
following Dirac equation [27, 29, 75, 77]: 

{'fd li + 7 5 (d z + 2d z A) - 7]e A F^) } * = 0, (4.2) 

with 7^9^ the 4-dimensional Dirac operator. 

It is convenient to decompose the full 5-dimensional spinor by 

*(X,Z) = ^2?pLn{x)fLn(z) + ^2^Rn{x)fRn(z). (4.3) 

n n 

Here iphn( x ) and ipRn(x) are the left-hand and right-hand components of a 4-dimensional 
Dirac field, which satisfy the four-dimensional Dirac equations 

l^d^nix) = m n ip Rn (x), 7 M 9 M V'Rn(x) = m n ijj Ln (x). (4.4) 
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Now, we are left with the following coupled equations 

[8 Z + n e A F(c/ ) )] f Ln (z) = m n f Rn (z), (4.5a) 
[d z - r, e A F(c/ ) )] f Rn {z) = -m n f Ln (z). (4.5b) 

Provided that eqs. (4.5a)-(4.5b), along with the following orthonormal conditions 

/oo 
e 4A fhm(z)fR n {z)dz = 5 RL 5 mn , (4.6) 
-oo 

are satisfied, we can recast action (4.1) into the standard four-dimensional actions for the 
massless and a series of massive chiral fermions. By redefining fi, n = e 2A /l„, we can 
rewrite eqs. (4.5) as two Schrodinger equations 

( - d 2 z + V L {z)) kn = m 2 J Ln , (4.7a) 
( - d 2 + V R (z))f Rn = m 2 J Rn . (4.7b) 

The effective potential for fi,n has been denoted as 

V L (y) = e 2A r] 2 F 2 (4>) - e A r)d z F((j>) - (d z A)e A V F(4>). (4.8) 

For f Rn , the coresponding potential V R (y) can be obtained from eq. (4.8) by simply chang- 
ing r] — > —i]. Note that in z-coordinate (defined by eq. (3.5)) 

4>{z) = ± - ==, 4.9 

^3 + 2X^)v 2 z 2 

A(z) = ~ In (l + ^X^v 2 z 2 ^j . (4.10) 

Thus if we want to obtain Z<i symmetric potentials Vl(?/) and V R (y), we have to constrain 
F((j)) to be an odd function of 4>. 

Now, our task is to solve the Schrodinger-like equation (4.7) with different forms of 
F(<j)) which determines the shapes of the potentials Vl and V R . In braneworld models we 
mainly interested in two types of effective potentials. 

The first one is the so called volcano potential, which has been discussed in section 
3 (see eq. (3.12)). For such potential we usually get a unique bound state, i.e., the zero 
mode, and a gapless continuum of KK modes. Sometimes, we might obtain several quasi- 
localized massive states (the resonant states). In section 3, we got no sign which indicates 
the existence of a resonant state. But as we will state below, such resonant states might 
exists for particular F(cj)). 

The other one is the well known Pdschl- Teller potential, which approaches to a positive 
constant rather than zero as z — > ±oo. This potential is very interesting, because in 
addition to the zero mode, there might be some massive KK states. Besides, the continuous 
KK modes are separated from the bound states by a mass gap. 

To proceed on, let us consider the following two cases, the first one leads to a volcano 
potential, while the other one a Poschl- Teller potential. Due to the fact that left- and 
right-handed fermions usually share similar resonant and massive spectrum [62], we only 
discuss the left-handed fermions. 
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4.1 Massive fermions as resonant states 

A we obtain a Schrodinger equation with the following potential 



By takeing F( 



2A< 5 ) 



3v 2 r](v 2 z 2 (3r] + 2X^) -3) 

(3 + 2v 2 z 2 A( 5 )) 2 
U7i] (U7i]z 2 + 29 (/cV - l) 4) 
841(1 + k 2 z 2 ) 2 4 



(4.11) 



which satisfies VL(0) = —v 2 rj and Vl(oo) = 0. That means only the left-hand fermion zero 
mode could be trapped on the brane for positive coupling constant r\. The zero mode, 
according to the discussions in [27, 77], takes the form 



dtuo(z) oc exp 



-V 



dz'e A ^'\ 



3 + 2A( 5 Vz 2 



(4.12) 



where (3 = |^5j is a dimensionless constant. It is easy to proof that ai,o(z) is normalizable 
if an only if 77 > 770 = A^ 5 ) /3. 

Since usually the fermions are massive, we hope to find a way to trap the massive 
fermions. For volcano potential, the massive bound state is absent. However, there might 
be some resonant states, which have large relative probabilities near the location of the 
brane. The existence of such resonant states offer us an effective way to trap massive 
fermions. This matter trapping mechanism is refereed to as the quasi-localization. Mimic 
to the case of gravity, we define 



^l,r(™ 2 ) 



$t h Zb \h,R{z)\ 2 dz 



nr i/l,r(^)i 2 ^' 

• J ^raax 



(4.13) 



as the relative probability for finding a left- (or right-) hand fermion resonant state fi,(z) 
(or /r(z)) with mass m. Sharp peaks in figure of i\,,R(m 2 ) can be identified as the resonant 
states. The resonant states are unstable, their lifetime can be estimated by r ~ T -1 , with 
r = 5m the width of the half height of the resonant peak. As shown in figure 5, for larger 
value of 77, there are more resonant states. The mass, width and lifetime of the resonant 
states are listed in table 1. 

4.2 Massive fermions as bound states 

A coupling given by 



tanh 



v 2 — <f) 2 

would lead us to the Poschl- Teller potential 

V L (z) = k 2 [(kr]) 2 - kr}{\ + kr])sech 2 (kz)] , 



(4.14) 



(4.15) 



- 11 - 



1 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 

- 



even 
odd 



10 20 30 40 50 60 70 80 



1 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 

- 



(i 



even 
odd 



50 



100 

™2 



J 50 



200 



Figure 5. The figures of Pl with k = 1, k§ = 1, r\ — 3 (the left panel) and rj = 5 (the right panel). 
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T 


3 


58.033367 


1 


27.28253 


5.22327 


0.000095 


10525.6787 






2 


47.73727 


6.90922 


0.0275038 


36.3586 






1 


47.62055 


6.90076 


2.43773 x 10~ 10 


4102176995.31678 


5 


160.820587 


2 


88.91114 


9.42927 


0.000005757 


173695.37532 






3 


123.32123 


11.10501 


0.0017165 


582.58201 






4 


149.47803 


12.22612 


0.0540928 


18.48673 



Table 1. The mass, width, and lifetime of resonances of the left-handed fermions (shown in 
figure 5). The parameters are k = 1, n 



1. V^ lax is the maximum of the height of the potential 



Vh, and n is the order of resonant states corresponding m from small to large. 



for which Vl(0) = —k 3 rj, and Vl(oo) = k 4 r] 2 is a positive constant. Again, the zero mode 
exists only when 7] > 0: 



«L0 



(z) = JV cosh-*"(A:z). 



(4.16) 



&lo( z ) is normalizable, provide that kr\ > 1/2. For simplicity, let us take krj = 3/2 1 , then 
the potential (4.15) reduces to 



Vl(z) 



9 



15 



/c 2 sech(/cz) 2 . 



As we have studied in [29, 



4 4 

such potential supports two bound states: 



o:lo(z) oc cosh 3//2 (kz) 
ahi(z) oc sinh(A;z) cosh 3 / 2 (£;z) 



with 
with 



m 
m 2 



0, 

2k 2 . 



(4.17) 



(4.18) 
(4.19) 



1 In fact, when krj — 1,2, ... is a positive integer, the quantum system possesses very peculiar property, 
namely, the system is reflectionless (and such property turns out to be intimately related to nonlinear 
integrable KdV system). Such peculiar property of Poschl- Teller system Vl is reflected in the presence of 
a hidden nonlinear bosonized supersymmetry of order 2krj + 1 in it, see [78, 79]. We would like to thank 
Mikhail S. Plyushchay for reminding us of their interesting works. 
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We got only one massive bound state mode, i.e., cxli(z). However, as the product krj 
increases, more and more massive bound states appear. The number of the bound states 
equals to the integer part of kr\ [26]. In addition to the bound states, there is also a 
continuum of KK states which start from m 2 = k 4r q 2 and asymptote to plane waves as 
z> 1. 

5 Conclusions and prospect 

In this paper, we study a thick domain wall solution in squared curvature gravity. Different 
from the Gauss-Bonnet gravity, the dynamical equation in our model contains higher order 
derivatives. As hoped, the solution itself is smooth and regular, and is stable under the 
gravitational tensor perturbations. The localization of gravity and spin- 1/2 fermions is 
shown to be possible for our solution. For gravity, only the zero mode can be localized. 
While for the left-handed fermions, in addition to the zero mode, some massive KK modes 
can also be localized as bound states (with infinite life), or be quasi- localized as resonant 
states (with finite life). The number of massive KK modes depends on the thickness of the 
brane k ~ 7 _1 , and the strength of the scalar-fermion coupling, i.e., rj. The KK spectrum 
of the right-handed fermions is similar to the one of the left-handed fermions, except the 
absent of the zero mode. 

A new feature of our model is that the curvature square term might influence the 
thick of the wall via the relation k ~ 7 ■ As the temperature effects is considered, the 
coefficient 7 might be a function of temperature T, i.e., 7 = j(T). Then, the thickness of 
the brane can be tuned by temperature. Some massive KK states of the fermions might 
appear or disappear, for k is another factor which influences the number of the KK states 
of the fermions. Besides, the Mexican hat potential (2.6) indicates that there might be a 
spontaneous symmetry breaking, which leads to phase transition and finally the formation 
of our four-dimensional domain wall. The quantitative calculation is a non-trivial task, we 
hope some related works can be reported in the near future. 
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